Chirality tunneling in mesoscopic antiferromagnetic domain walls 
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We consider a domain wall in the mesoscopic quasi-one-dimensional sample (wire or stripe) of weakly 
anisotropic two-sublattice antiferromagnet, and estimate the probability of tunneling between two 
domain wall states with different chirality. Topological effects forbid tunneling for the systems 
with half-integer spin S of magnetic atoms which consist of odd number of chains n±. External 
magnetic field yields an additional contribution to the Berry phase, resulting in the appearance of 
two different tunnel splittings in any experimental setup involving a mixture of odd and even n±, 
and in oscillating field dependence of the tunneling rate with the period proportional to l/n±. 
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In recent years, there have been much interest to the 
problem of quantum spin coherence in mesoscopic mag- 
netic systems, mainly in nanoparticles Experimen- 
tally, one of the crucial problems is to prepare an en- 
semble of particles of highly uniform size: even small 
fluctuations of size lead to large fluctuations of the tun- 
neling probability since they contribute to the power of 
exponent. One nice solution is to use high-spin molecu- 
lar clusters another possible way, proposed in PLQ, 
is to use topologically nontrivial magnetic structures: do- 
main walls in quasi-one-dimensional (Id) systems (wires, 
stripes), vortices in 2d systems, etc. Such objects have 
mesoscopic scale, e.g., in weakly anisotropic magnets 
with magnetic ions in s-states the domain wall thickness 
is usually about 100 lattice constants, and since their 
shape is determined by the material constants they are 
to a high extent identical. 

Classically, magnetic domain wall (DW) has certain 
"chirality," an internal degree of freedom characterizing 
the way of rotation of magnetization inside a DW. Two 
states with opposite chirality are completely equivalent 
in energy (we will not consider magnets without inversion 
center where this is not true). In quantum case there is 
generally a nonzero transition amplitude mixing the two 
states and lifting the degeneracy j^Pl; that under favor- 
able circumstances this tunnel splitting can be detected 
with resonant technique of some kind. It is known that in 
antiferromagnets (AF) tunneling is generally more favor- 
able than in ferromagnets, both in case of fine particles 
H and domain walls B . 

In this paper we show that in the simplest model of 
mesoscopic AF with half-integer spin S of magnetic ions 
topological effects forbid chirality tunneling for a DW 
with odd number n± of spins in its cross-section. Wc 
further show that in presence of even weak external mag- 
netic field this strict "selection rule" is relaxed, which 
leads to appearance of two different values of tunnel split- 
ting in any half-integer S sample with weakly fluctuating 
nj_. For any S, the tunneling amplitude is shown to be 
an oscillating function of the field. 



Model. Consider a thin quasi-one-dimensional stripe of 
two-sublattice weakly anisotropic antiferromagnet, which 
we for the sake of simplicity consider as a system of n± 
AF chains of spin-S" magnetic atoms, coupled antiferro- 
magnetically with the same exchange constant J for any 
neighboring spins. We assume that magnetic atoms form 
a perfect crystal structure on a bipartite lattice, as shown 
in Fig. |l|; note that nj_ can be odd or even without intro- 
ducing any frustration. We assume a rhombic anisotropy 
of the form 



(1) 



where i labels lattice sites and K\p "C J are the 
anisotropy constants, K\ > K 2 > 0, so that OZ is the 
hard axis and OX is the easy axis in the easy plane XY. 

Due to the quasi-ld structure, one can assume that 
the sublattice magnetization depends only on the space 
coordinate x along the wire (note that X and x axes do 
not need to coincide). Using the standard technique ||, 
one can obtain the effective Euclidean action of AF in 
continuum approximation, which has the form of a well- 
known 0(3) nonlinear sigma-model: 

A E = ~n±hSW + i2<!m±Sh(Q + Q' H ), (2) 

w[i] = J d 2 x[{dj){dj) + ^[{i + P )i 2 z + i 2 Y } +w a }, 

Q= — [ d 2 xl- (d x l x d 2 l) , 
4tt J 

Q'h = ^- c J d 2 x H -(lxd 2 l). 

Here I is the unit Neel vector, {x\,X2) — {x,ct) is the 
Euclidean plane, c = JSaZ c /h is the limiting velocity 
of spin waves, Z c is the lattice coordination number, a 
is the lattice constant, A = a(J Z c / AK2) 1 / 2 ^ a is the 
characteristic DW thickness, p — KxjK% — 1 is the rhom- 
bicity parameter, 7 = gpg/h is the gyromagnetic ratio, g 
denotes the Lande factor, and [1b is the Bohr magneton. 
The quantity w a {l) — (j/c) 2 (H ■ I) 2 describes effective 
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renormalization of the anisotropy induced by the field. 
In (j|), the term proportional to Q is the so-called topo- 
logical term originating from the sum of Berry phases 
of individual spins, Q being the homotopical index of 
mapping of the (cci , x-i) plane onto the sphere I 2 = 1, and 
Q' H is the contribution from magnetic field. 

A static DW solution Iq(x) corresponds to the rotation 
of vector I in the easy plane XY: 

lox = c'tanh(a;/A), l y = cosh(x/A), l oz = , (3) 

where a, a' — ±1. The quantity a' is the "topological 
charge" of the DW, and the chirality a determines the 
sign of I projection onto the "intermediate" axis OY . 
Two states with a = ±1 are equivalent in energy; change 
of a describes reorientation of the macroscopic number 
of spins N DW - A/a > 1, typically N DW ~ 70 4- 100. 

Chirality tunneling in absence of magnetic field. Let 
us consider first the case H = 0. Tunneling between the 
DW states with opposite chiralities can be studied using 
the instanton formalism. Since the tunneling here occurs 
between two inhomogeneous states, the corresponding 
instantons are non-one-dimensional (space-time). The 
structure of instanton solution li us t{x, r) is shown in Fig. 
^; it has the following asymptotic behavior: 

lx — > ±c', x — > ±oo, ly — > =pcr, x — 0, t — » ±oo 
J z ->p = ±l, a?->0, r->0, (4) 

note the appearance of another topological charge p = 
±1. Along any closed path in the Euclidean plane go- 
ing around (but far from) the instanton center vector I 
rotates by the angle 2ttv in the easy plane XY, where 
v = era' = ±1. Thus, the instanton configuration has 
the properties of an out-of-plane magnetic vortex (i.e., 
with lz in the center) and is characterized by two 
topological charges ||: vorticity v and polarization p. 
The instanton solution satisfies the equations 

V 2 (9 + sin(9cos6l[(p + cos 2 ip)/A 2 - (V<^) 2 ] = 0, 

V • (sin 2 OVif) - (1/ A 2 ) sin 2 sin ip cos ip = 0, (5) 

where we have introduced angular variables as Ix+Hy — 
s'mOe^, lz — cos 6*, and V = (81,82) denotes the Eu- 
clidean gradient. 

We are not able to find the exact solution of Eqs. (||), 
but the tunneling amplitude can be estimated from ap- 
proximate arguments. One can readily observe that, in 
contrast to the same problem for ferromagnet ||, Eqs. 

, as well as their solutions, are real. Thus in absence 
of magnetic field the real part of Ae is given by n±SW/4, 
and the imaginary part is completely determined by the 
topological term Q (note that Q is a total derivative and 
does not contribute to the equations of motion). Then 
the procedure of constructing the instanton solution can 
be viewed as minimization of W; it should be remarked 



that the real part of the Euclidean action for the instan- 
ton formally coincides with the energy of a vertical Bloch 
line in a 2d DW, which is rather well studied [O. 

Another observation is that in AF in absence of mag- 
netic field the translational DW motion in real space and 
its internal degree of freedom (chirality) are completely 
uncoupled, in contrast to the situation in ferromagnets 
H ; for H ^ this coupling appears but becomes impor- 
tant only in strong field regime . Thus for weak fields 
one can calculate the transition amplitude (^|) assuming 
that the DW coordinate in real space is just fixed. 

Further, for uniform boundary conditions at infinity 
the quantity Q determining Im.4^; can take only integer 
values, but in our case Q = —pv/2 = ±^ is half-integer, 
which is typical for out-of-plane vortices (see, e.g., ||). 
For the given a 1 there are two instanton solutions with 
the same vorticity v and opposite polarizations p which 
equally contribute to Re.4E but have different signs of 
\yclAe- The tunnel splitting is given by 

T = CfuJi { nA WS/^) 1/2 e~ ni -W s/i \ cos$| , (6) 

where W = W[l- mst {x, r)] — Fy[Z (a;)] is the difference 
of the real part of Euclidean action calculated on the 
instanton solution and on a static DW without instan- 
ton, u>i = (c/A)y/p is the frequency of the out-of-plane 
magnon localized at the DW playing the role of "attempt 
frequency," C is a numerical constant, and $ is the phase 
determined by the p-dependent part of Im^^. In the 
simplest model with H = considered so far $ = irn±S, 
and thus the tunneling amplitude vanishes when S is half- 
integer and n± is odd. In any mesoscopic sample with 
weakly fluctuating cross-section the value of n± » 1 will 
be randomly odd or even for different samples or different 
parts of the same sample. Thus the tunneling is forbidden 
for approximately one half of domain walls in half-integer 
S antiferromagnets. This is somewhat similar to the case 
of half-integer S nanoparticles with non-compensated to- 
tal spin S*tot where the Kramers theorem forbids tunnel- 
ing for the particles with half-integer Stot [@ ■ Below we 
will see that in presence of external magnetic field the 
situation is different: both integer and half-integer n± 
contribute, but with different tunneling rates. 

One can observe that the problem has three different 
length scales: the DW thickness A, the vortex core size 
in spatial direction A ra , and the characteristic size of 
the core in the imaginary time direction A VT . For strong 
easy-plane type anisotropy, p ^> 1, the vortex core is 
nearly axially symmetric: up to distances r <C A the 
anisotropy in the easy plane can be neglected, and the 
solution in the core reduces to the well-known case of a 
usual vortex in an easy-plane magnet Jl3[ , with 8 — 6q(£), 
<P = vx, v = ±1, 

d 2 6 /d£ 2 + (1 - v 2 H 2 ) sin O cos Q = , (7) 

here (r, x) are the polar coordinates in (xi,X2) plane, 
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r = [x\ + x 2 ) 1 / 2 , % = arctan^/xi), and £ = r^fp/A. 
Thus one has A vx = A VT = A/^/p <c A, so the core is 
isotropic and much smaller than the DW thickness. 

In the opposite "almost easy-axis" case p<C 1 the vor- 
tex core is strongly asymmetric: its spatial size A vx coin- 
cides with the DW thickness A, but the imaginary time 
size A VT = A/y/p is much larger 

On the other hand, for r 3> A VXl A VT , i.e., far outside 
the core, one can put 9 ~ 7r/2, which reduces the system 
(|) to the 2d elliptic sine-Gordon equation, 

VV = (l/2A 2 )sin2^. (8) 

In the large rhombicity limit p 1 within a wide range 
of r (for A/^fp < r < A) the solutions (@) and (§) 

can be regarded as coinciding, and the integrand in W 
is proportional to 1/r 2 . Then, one may divide the in- 
tegration domain into two parts: r < R and r > R, 
where R is an arbitrary point in between A/^fp and 
A. For r < R the solution (]?]) may be used, yielding 
W r< R = 2ir\n((Ry/p/ A) with the numerical factor £ be- 
ing known from Ref. |l3), ( ~ 4.2. For r > R, one can 
use a trial function approximately satisfying (®), e.g., 

tan</3 = tanh[a; 2 /(A~ ;I . 1 A„ T A)]/ sinh(xi/A) , 

cosO = [cosh(a;i/A ra ) cosh^/A,,,-)] -1 , (9) 

and evaluate the integral in W numerically, which for 
p > 1 gives W r>R = 27rln(C'A/i?) with (' ~ 0.525. 
Summing up the two contributions, we obtain 

W ~ 27rln(2.2Vp), P>1. (10) 

In the weak rhombicity limit p <C 1 the trial function 
(^) can be used for the entire (x\,X2) plane, yielding the 
result 

W~8p 1/2 , p<l, (11) 

which coincides with one obtained by us earlier in the 
effective Lagrangian approach The result ( |ll| ) 

breaks down only for extremely small rhombicities p <C 
4/(n±S) 2 , when the system is very close to the easy- 
axis regime and its low-energy spectrum coincides with 
that of the free rotator, which yields the tunnel splitting 
r ~ {hc/2An ± S) fL§. 

Comparing the results for tunneling in AF domain wall 
with those for a ferromagnetic DW M, one can see that 
(i) for a ferromagnetic DW the function W determining 
the tunneling exponent contains_the additional large fac- 
tor A/a; (ii) for ferromagnet W oc J~p at large p while 

for AF the growth of W(p) is much slower. 

Effect of magnetic field. Consider now the behavior of 
the imaginary part of the Euclidean action when a weak 
external magnetic field H is applied to the system (we 
neglect here the anisotropy renormalization w a because 



its effect is rather trivial). One can see that the mixed 
product in Q' H significantly differs from zero only in the 
vortex core, and thus in case of large rhombicity p ^> 1 
the isotropic vortex solution (0) may be used to estimate 
it. After integration we obtain 

Q' H ^ P X(H X /H C ), (12) 

where H c = {AZ C S 2 JKf) 1 ' 2 /-f denotes the magnitude 
of field for which the field-induced anisotropy becomes 
equal to the easy-plane one, Hx is the field component 
along the easy axis, A = J Q d£{isin2# + £(d6o/d£)} 
is a numerical constant, A ~ 3.83. This results in the 
following expression for the phase factor $ in (||): 

$^$ H =im ± S[l + (\H x /2H c ], (13) 

thus the tunneling amplitude oscillates as a function of 
the field with the orientation-dependent period 

SH = (2H c /n ± SX) . (14) 

This period can be rather small: assuming 5 = 1 and a 
typical H c - 100 kOe, one gets 8H ~ 2 kOe-=-20~ Oe for 
n± = 10-i- 10 3 . A similar oscillating behavior was pre- 
dicted earlier for tunneling in small ferromagnetic |l6| ] 
and antiferromagnetic Jl7],|l8| particles, with the differ- 
ence that in the AF case instead of the field H c tx V 'JK\ 
in (Jl4|) a much stronger exchange field H e oc J would be 
present. For half-integer S presence of the magnetic field 
lifts the degeneracy of two odd-?ij^ DW states with op- 
posite chirality, allowing tunneling between them. Note 
that the "correction" in $ is strong even for weak H since 
it contains a large factor n±. 

Another consequence of the above result is that for 
half-integer S there are two different values of the tun- 
nel splitting for even and odd n± , which means that for 
any experimental setup with weakly fluctuating n± there 
should be two different resonance peaks at the corre- 
sponding frequencies; this effect was overlooked in previ- 
ous studies. It is worthwhile to remark that the same ef- 
fect should be also present in half-integer S AF nanopar- 
ticles with noncompensated spins considered first by Loss 
et al. jL2| and also studied later in Refs. M], provided 
that there is some weak interaction (e.g., magnetic field 
or the Dzyaloshinskii-Moriya interaction |p0[ ) contribut- 
ing to the phase factor $ and shifting it from a multiple 
off. 

We would like to finish with a word of caution: in pres- 
ence of field the problem of chirality tunneling is actually 
more complicated then one can guess from the simple ar- 
guments presented above. The point is that the field 
contribution Q' H , unlike Q, is not a total derivative and 
thus yields an imaginary perturbation to the equations 
of motion, causing nontrivial changes in the instanton 
structure and in the spectrum of eigenmodes in presence 
of instanton which eventually contribute to the phase fac- 
tor $ JTq]. However, using a perturbation theory in H, 
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one can show that corrections from the change of instan- 
ton structure contribute to W as (H/H c ) 2 and to <f> as 
n±(H I H c ) 3 , so that they can be neglected for weak fields, 
and the contribution to $ from fluctuations does not con- 
tain n±. Thus we think that the formula ( |l4| ) remains 
correct, together with our main conclusions on the os- 
cillating field dependence of T and on presence of two 
different tunneling rates. 
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FIG. 1. (a) a schematic picture of the cross-section of an- 
tiferromagnetic mesoscopic stripe; (b) two domain walls with 
opposite chiralities. 
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FIG. 2. The structure of instanton connecting two DW 
states with opposite chiralities. Arrows denote projections of 
vector I on the easy plane, and on the thin solid line vector I 
forms the angle of about 45° with the easy axis. 
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